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We evaluate the ttN coupling constant using a QCD sum rule based on the 



I/") ■ pion-to-vacuum matrix element of the correlator of two interpolating nucleon 

ON 

r*| ■ fields. The part of the correlator with Dirac structure ^75 is used, keeping all 



terms up to dimension 5 in the OPE and including continuum contributions 



■ °n the phenomenological side. The ratio of this sum rule to the nucleon sum 



rule involving condensates of odd dimension yields stable results with values 



of g-wN in the range 12 ± 5. The sources of uncertainty are discussed. 

Understanding hadron interactions from first principles is one of the main goals of Quan- 
tum Chromodynamics (QCD). Since the solution of QCD for hadron interactions at low 
energies is still far off it is useful to consider tackle problems of hadron dynamics with ap- 
proaches that lie as close as possible to QCD. One of them, the method of QCD sum rules 
H, has proved to be a very powerful tool to extract information about hadron properties. 
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We present here a sum rule analysis of the ttN coupling constant g-^N, which is one of the 
most important parameters in hadron physics. This quantity has previously been examined 
within the framework of QCD sum rules in Refs. Reinders, Rubinstein and Yazaki 

worked with two different sum rules for g n N, one obtained from a correlator of three 
interpolating fields, and one based on the pion-to-vacuum matrix element of the correlator 
of two interpolating nucleon fields, 77: 

(0\T{r,(x)fj(0)}K(k)), (1) 

However those authors included only the leading term of operator product expansion (OPE) 
and neglected continuum contributions. Shiomi and Hatsuda developed the sum rule 
based on the part of this two-point correlation function with Dirac structure 75. They 
calculated g n N in the soft-pion limit, including condensates up to dimension 7 in the OPE, 
radiative corrections and a perturbative estimate of continuum contributions. 

The method based on the two-point correlator ([I]) has a significant advantage in that 
it can be used for low values of the momentum transfer to the nucleon. In contrast the 
OPE of the three-point correlator is valid only for large spacelike meson momenta and so a 
determination of the coupling constant requires an extrapolation to zero momentum (see, for 
example [|J). This procedure is dangerous because of the higher-order terms that have been 
omitted from the OPE. The contributions of these terms give rise to corrections that are 
proportional to large inverse powers of the meson momentum k, making the extrapolation 
of a truncated OPE unreliable. Estimates of the coupling constant from the coefficient of 
1/k 2 determined at large k 2 , as in Refs. |5|||, cannot distinguish the meson pole term from 
the contributions of higher-mass states in the same meson channel. 

In the soft-pion limit the OPE for the 75 part of the two-point correlator for g^N has 
exactly the same form as that for the nucleon sum rule |||7| involving condensates of odd 
dimension, up a factor of 1/ f w [@|§. If continuum corrections are neglected, the ratio of these 
two sum rules has the form of the Goldberger-Treiman relation with g^ = 1 • Shiomi and 
Hatsuda |§ showed that this continues to hold for the higher- dimension terms in the OPE, 
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provided that the continuum thresholds are taken to be the same for both sum rules. Using 
different thresholds in the two sum rules they were able to get around this problem with the 
implied value of gA- 

However, using the usual soft-pion theorem ||, the correlator (1) can be expressed as 

"f (0|[g^T(r ? (x),r ? t(0))]|0) = -^{ 75 r a ,(0|T( ?7 (x),r ? t(0))|0)} (2) 

where Q% is the axial charge and we have made use of the transformation properties of 
the interpolating field under axial rotations ||, [Qf,^] = — |7sT a ^- The anticommutator 
with 75 picks out the part of the two-point correlator proportional to the unit Dirac matrix. 
The phenomenological side of the resulting sum rule is thus 275 / f m times the corresponding 
expression for the odd-condensate nucleon sum rule. This matches exactly with the structure 
found for the OPE side in Refs. gfj. 

The soft-pion limit for the 75 piece of the correlator (1) thus yields a sum rule for 
Mn/ f n = Q-kn/qa- The value for g nN determined from such a rum rule follows from the 
odd-condensate sum rule for the nucleon mass and the Goldberger-Treiman relation (or an 
approximation to it taking qa — 1). The sum rule can be thought of as just a chiral rotation 
of the odd-condensate nucleon sum rule and not an independent determination of g n N. By 
considering terms beyond the soft-pion limit, we obtain here a value for g^N that is not 
simply a consequence of chiral symmetry. The sum rule we use is thus analogous to that for 
the 7riVA coupling || and the approach can be applied to calculations of other meson-baryon 
couplings. 

Moreover, a potentially important piece of the phenomenological side of the sum rule for 
g n N has been omitted in both calculations. This term corresponds to transitions of where 
a ground-state nucleon absorbs the pion and is excited into the continuum. Since it is not 
suppressed by the Borel transformation such a term should be included in a consistent sum- 
rule analysis, as pointed out long ago |7|{n| and stressed recently by Ioffe [[TI|,[nj. In the 
soft-pion limit, such terms generate contact interactions where the pion couples directly to 
the nucleon field, (N(p)\TJ n (0)\Tc(k)), and which are essential if the correct soft-pion limit 
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is to be obtained. Although the need for these terms is particularly clear if pseudovector 
ttN coupling is used (see for example |I3"|), they should also be included for pseudoscalar 



coupling. Their omission in Refs. ||[J can explain why the correct Goldberger-Treiman 
relation was not found there. Indeed, as the authors of || point out, a quick estimate of 
these unsuppressed N* contributions suggests that they could be as large as 25%: enough 
to remove the discrepancy with the Goldberger-Treiman relation. 

Here we start from the two-point correlator ([]]) just discussed, but instead of the piece 
with with Dirac structure 75 considered in Refs. |||| we work with the structure ^75, where 
k is the pion momentum. We work here to leading order in a chiral expansion, neglecting 
higher-order terms in the pion momentum or current quark mass. 

We consider the two-point correlation function 

n(p) = i J d 4 xexp(ip ■ x)(0\T{ Vp (x)rj n (0)}\7r + (k)), (3) 
where we use the Ioffe interpolating field || for the proton, 

V P (x) = e afec [^(x) T C7X(x)]757^ c (^, (4) 

where a, b, c are the colour indices and C is the charge conjugation matrix. (The corre- 
sponding neutron field i] n is obtained by interchanging u and d quark fields.) More general 
choices of interpolating field are possible, as discussed in detail by Leinweber |L4j]. For the 



odd-condensate nucleon sum rule, which we make use of in our determination of g n N, it 



turns out that the Ioffe field is close to the optimal one as determined in Ref. JTj] and so we 
do not consider more general fields here. 

In the deeply Euclidean region {p 2 large and negative) the OPE of the product of two 
interpolating fields takes the following general form 

i U x ^ii P . x )T { r, p ( x )nM} = ^CMO„ (5) 

J n 

where C n (p) are the Wilson coefficients and O n are local operators constructed out of quark 
and gluon fields (all renormalised at some scale fx). Using this OPE in the correlator (|j), we 



find that only operators of odd dimension contribute. The leading term in this expansion 
involves operators with dimension 3 and has the form 

n 3 (p,fc) = -^p 2 ln(-p 2 )(0|S 7 Q 75M|7r + (A ; ))7 a 7 5 + ---, (6) 

where terms that do not contribute to the Dirac structure of interest, ^75, have been sup- 
pressed. The matrix element here is just the usual one for pion decay: 

(0\d 1 a l5 u\n + (k))=iV2f 7T k a , (7) 

where f n = 93 MeV is the pion decay constant. 

At dimension 5 the only relevant contribution arises from the second-order term in the 
covariant expansion of the nonlocal operator <i(0) 7 a 7 5-u(x). It has the form 

n 5 (p) = ^H-p 2 )(0\d 1 a l5 D 2 u\n + (k)) 7al5 + ---. (8) 

Up to corrections of higher order in the current mass, this can easily be re-expressed in 
terms of a mixed quark-gluon condensate 

(0|3 7 Q 75J D 2 M |7r+(£;)> = ^-(0\d 7 a 7 5 cr^G^u\n + (k)) + 0{m\). (9) 

Some further manipulation allows one to rewrite this in the form 

(0|3 7 Q 75J D 2 u|7r+(£;)> = -g s ((Q\dG a ^u\7r + (k)) - ig s (Q\dG^ 5 u\7r + (k))), (10) 

where G^ v = ^e^ upa G pa . (We use the convention e 0123 = +1.) The second term in this 



expression is of higher order in the chiral expansion (see Ref. [15] for details) and so we 
neglect it here. 

The first matrix element in (|10[) was extracted by Novikov et al. |15j from two QCD sum 
rules for the pion. They expressed it in the form 

g s (0\dG a ^u\n + (k)) = V2i6 2 f n k a , (11) 
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and obtained 6 2 = (0.20 ± 0.02) GeV 2 .[] A crucial contribution in both of their sum rules 
is the four-quark condensate, a s (0\(qq) 2 \0). Novikov et al. |15j used the factorisation ap- 



proximation for this quantity but direct determinations of it from other sum rules suggest 
significantly larger values |T6|-|18||, at least 2-3 times bigger than those obtained from fac- 
torisation. These give correspondingly larger values for 5 2 , a point we shall come back to in 
the discussion of our results below. 

As an estimate of the importance of higher dimension condensates, we have also calcu- 
lated the contribution of what we hope is the most important dimension-7 operator in the 
OPE. This is a mixed quark-gluon condensate, which we evaluate in the factorised approx- 
imation. Keeping only this contribution explicitly, the dimension-7 piece of the correlator 
is 

1 _ n 
n 7 (p) = - T ^(0|ci 7 a 75 n|7r + (fc))(0|-iG 2 |0) 7Q 75 + ■ ■ ■ , (12) 

where (0|^G 2 |0) is the gluon condensate in vacuum. 

On the phenomenological side, the term of interest in the correlator @ is the one with 
a double pole at the nucleon mass, since this contains the ttN coupling constant. However 
the nucleon interpolating field does not just create ground-state nucleons; there are also 
continuum contributions which cannot be ignored. The continuum-to-continuum pieces are 
modelled in the usual manner, in terms of the spectral density associated with the imaginary 
part of the OPE expression for the correlator. This continuum is assumed to start at some 
threshold SVtv- After Borel transformation, it can be taken over to the OPE side of the sum 
rule where it modifies the coefficients of the terms involving ln(— p 2 ). In addition one must 



x There is a potential sign ambiguity in using the result of Ref. [15] since they do not specify 
their convention for e 0123 . However we have checked that our coefficient of x 2 in the expansion of 
(0|d(0)7 a 75ti(x)|7r + ) (which should be independent of convention) agrees in both sign and mag- 
nitude with that of the corresponding term in the expansion of (0\u(x)u(0)\0) a given in Refs. 
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include nucleon-to-continuum terms since Borel transformation does not suppress these with 
respect to the double-pole term p|,|IO|-|T2| . To first order in k, the correlator has the form 



n(p) = iV2# 



75 



X 2 N M N g wN r°° 1 / 1 a(s) 



dsb(s)- —r [- -=- + 



+ •••, (13) 



> 2 - M 2 N ) 2 Jw* v ' s- M 2 N \p 2 - M% s-p 2 1 

where the continuum-continuum terms (and terms with other Dirac structures) have not 
been written out. Here A at is the strength with which the interpolating field couples to the 
nucleon: 

(0\ m (0)\N{p)} = \ N u(p). (14) 

Note that the strength of the ^75 piece of the double-pole term is the same for both pseu- 
doscalar and pseudovector ttN coupling. 

Equating the OPE and phenomenological expressions for the correlator (Q) and Borel 
transforming [||], we get the sum rule 

1 r M 4 £ 2 (x) + ^-M 2 E x {x)5 2 + ±-E (x)(0\^G 2 \0) = ( ^%f^ + A ) exp(-M^/M 2 ), 



2tt 2 zv ' 9vr 2 1V ' 12 v /N 1 7T 1 ' V f w M 2 



(15) 



where M is the Borel mass and E n (x) = 1 — (1 + x + ... + ^r)e~ x with x = ^f-. The 
second term on the r.h.s. of this sum rule, involving the undetermined constant A, is the 
Borel transform of the nucleon pole term of the nucleon-to-continuum piece in (|T3|) . It 
contains the same exponential as the nucleon double-pole term and so cannot be ignored. 
The second nucleon-to-continuum term in fllTf ) leads to a term that is suppressed by an 
exponential involving the masses of states in the continuum. It is thus typically a factor of 
3-4 smaller than the term included in (fL5|). Provided that the first of these mixed terms 
is a reasonably small correction to the sum rule, it should be safe to neglect the second, as 



discussed by Ioffe [111 



We now turn to the numerical analysis of this sum rule. First, one should get rid 
of the unknown constant A. Multiplying the sum rule by M 2 exp M^/M 2 , we see that 
the r.h.s. becomes a linear function of M 2 . By acting on this form of the sum rule with 



(1 — M 2 d/dM 2 ) |7] (or equivalently by fitting a straight line to the l.h.s. and extrapolating 
to M 2 = ||10|| ) we in principle can determine the value of gvjv- However we are unable to 
find a region of Borel mass in which the l.h.s. is approximately a linear function of M 2 , and 
hence there is no region of stability for the extracted g^N- 

This lack of stability is similar to the situation for the nucleon sum rules, where two 
sum rules can be derived || (involving either odd or even dimension operators) but neither 
shows good stability. Nonetheless the ratio of these leads to a more stable expression for the 
nucleon mass. We have therefore taken the ratio of our sum rule ( |i~5"D to the nucleon sum 
rules. We obtain the most stable results from the ratio to the odd-dimension sum rule, 

1 pM^^KOl^lO) + ^(0|gg|0)(0|^G 2 |0) = X 2 N M N exp(-M 2 N /M 2 ), (16) 



47r 2- xv-.wx-!^,-/ ■ 24 x-,«,-, x-, n 

and so we present here only the results for that case. Taking such a ratio also has the 
advantage of eliminating the experimentally undetermined strength A at from the sum rules. 
Note that we have allowed for a different continuum threshold Sn in the nucleon sum rule 
and have defined x^ = S^/M 2 . 

The ratio of the sum rules flT5| ) and ( |T6"D can be written in the form 

&M*E 2 {x) + £iM± El (x)S 2 + ±M 2 E (x)(0\^G 2 \0) 
U - i ^M^ 1 (x iV )(0|gg|0) + ^(0fe|0)<0|fG 2 |0) 9«n + A M , [U ) 

and the method discussed above used to eliminate the unknown mixed nucleon-to-continuum 
term, A'M 2 {A' = Af n / '\ 2 n Mn)- In Fig. 1 we show results for g^N as a function of M 2 , for 
typical values of the condensates and thresholds: (0|gg|0) = —(0.245 GeV) 3 , (0|^G 2 |0) ~ 
0.012 GeV 4 , 5 2 = 0.35 GeV 2 , S N = 2.5 GeV and S wN = 2.15 GeV. Stable values of g vN ~ 
11.7 are found over a region M 2 ~ 0.8 — 1.8 GeV 2 . Corrections due to the A'M 2 term are 
small, at most 5%. The second such term in (|T3| ) is expected to be smaller by a factor of 
3-4, and so we are justified in neglecting it. 

The threshold SVtv has been adjusted to give stability for M 2 around 1 GeV 2 , since one 
may hope that in this region the Borel transformed sum rule is not too sensitive to the 
approximations that have been made on both the OPE and phenomenological sides of the 
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sum rule. The existence of a window of stability provides a check on the consistency of this 
assumption. We also demand that the thresholds Sn and S^n should lie significantly above 
this window so that the continuum is not too heavily weighted in the Borel transform. We 
find that the window of stability moves rapidly upwards as S w n is increased for fixed Sn- 
For the typical parameter values above, only the region 2.05 GeV 2 < S^n < 2.22 GeV 2 
satisfies these requirements. The value of g n N varies by at most ±0.2 over this region. 

As a further check on our results, we have examined whether the individual sum rules 
(|15|) and fljj] ) satisfy the criteria suggested by Leinweber fly]]. We find that the highest 



dimension condensates contribute less that 10% of the OPE to both sum rules for M 2 > 0.8 
GeV 2 . The continuum forms about 40% of the phenomenological side of the differentiated 
version of the g n N sum rule ( |T5"D for M 2 up to 1.4 GeV 2 , the point at which the continuum 
reaches 50% of the odd-condensate sum rule (0). The region M 2 ~ 0.8 — 1.4 GeV 2 thus 
provides a window within which our results are both stable with respect to the Borel mass 
and not too sensitive to our approximations. 

We have examined the dependence of our results on the other input parameters. Variation 
of the threshold in the nucleon sum rule Sn from 2.2 to 2.8 GeV 2 , readjusting S n N to 
maintain stability, changes g n N by ±0.2. To estimate the sensitivity of our sum rule to the 
contributions of dimension-7 condensates and to uncertainties in the gluon condensate, we 
have varied the dimension-7 term in ( |i~5|) between zero and twice its standard value. Our 
results for g n N change by ±0.5 over this range. 

One of the most important input parameters in our sum rule is the matrix element <5 2 , 
defined by (pTTj) . As already mentioned, this parameter was extracted by Novikov et al. 
15| from an analysis of two sum rules for the pion. Their results depend crucially on the 



four-quark condensate, ct s (0|(gg) 2 |0), for which they made the factorisation approximation 
and took a value of about 2 x 10~ 4 GeV 6 . With this input, both of their sum rules yield 
consistent results for 5 2 in the region 0.20 ± 0.02 GeV 2 . However, sum rules analyses of 
r decay and e + e~ annihilation into hadrons lead to significantly larger values of the four- 



quark condensate (see fl6HT8]1 and references therein), in the range (4 — 6) x 10 GeV 



Using these in the sum rules of Ref. |15| leads to values for S 2 ranging from 0.28 to 0.45, 
although the two sum rules do not then give consistent results. As a conservative estimate 
of the uncertainty in 5 2 we have considered the range 0.20 to 0.45 GeV 2 . The corresponding 
variation in g n N is ±2, when the other parameters are held at their values above and S^n is 
changed to keep the window of stability around 1 GeV 2 . 

A second significant source of uncertainty is the quark condensate (0|gg|0) which appears 
in the odd-dimension sum rule for the nucleon. "Standard" values for this lie in the range 
-(0.21 GeV) 3 and -(0.26 GeV) 3 . The values of the nucleon mass determined from sum 
rules |J are strongly correlated with this condensate. There is also a weaker correlation 
with the chosen value of the threshold Sjy. Since we are dividing by M N in the ratio (|17D, 
our results are rather sensitive to the value of this condensate. One would like to use values 
of (0|gg|0) and Sn that give the nucleon mass correctly, but the ratio of the odd and even 
dimension nucleon sum rules does not yield completely stable results for Mjv- The best we 
can do is to rule out values of — (0|gg|0) below (0.23 GeV) 3 since they cannot reproduce the 
nucleon mass within the region of Borel mass and threshold that we consider. Varying the 
quark condensate between —(0.23 GeV) 3 and —(0.26 GeV) 3 , we find that g^N changes by 
±2. 

Our final result for the ttN coupling constant is thus g n ^ — 12 ±5, where the uncertainty 
is dominated by 5 2 and (0|gg|0). This is to be compared with values deduced from NN and 



ttN scattering. For many years the accepted value was g n N — 13.4 |19] but this coupling 



has been the subject of some debate in recent years. More recent analyses lead to values in 
the range 12.7-13.6 ||20|| . Our result is obviously consistent with any of these. The rather 
large uncertainty in it could be reduced if the quark condensate could be determined more 
precisely. In addition, the sum rules of Novikov et al. [jnj should be re-examined using larger 
values of the four-quark condensate to try to pin down S 2 more exactly. 

In summary, we have calculated the ttN coupling constant using a QCD sum rule based on 
the pion-to-vacuum matrix element of a two-point correlator of interpolating nucleon fields. 
This approach avoids the need for extrapolation from large spacelike meson momenta. We 
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have included nucleon-to-continuum terms omitted from previous analyses. Our sum rule 
is based on the part of the correlator with Dirac structure ^75 and includes all terms up 
to dimension 5 in the OPE. Stable results are obtained from the ratio of this sum rule 
to one for the nucleon mass and the unsuppressed nucleon-to-continuum contributions are 
found to be small. Contributions from higher-dimension operators and omitted continuum 
contributions are estimated to be small. This demonstrates the practicability of this type of 
sum rule for calculation of other meson-baryon couplings, whose values are at present not 
well determined. 
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FIGURE CAPTION 



Fig. 1. Dependence on the square of the Borel mass of the ixN coupling constant determined 
from the ratio of sum rules fllTD- The values of the parameters used are given in the text. 
The solid line shows the value of g n w corrected for the mixed continuum term A'M 2 , the 
dashed line the uncorrected l.h.s. of (|iT|). 
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